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The ANTARES telescope is well suited to detect neutrinos produced in astrophysical transient 
sources as it can observe a full hemisphere of the sky at all times with a high duty cycle. The 
background and point source sensitivity can be drastically reduced by selecting a narrow time 
window around the assumed neutrino production period. Radio-loud active galactic nuclei with their 
jets pointing almost directly towards the observer, the so-called blazars, are particularly attractive 
potential neutrino point sources, since they are among the most likely sources of the observed ultra 
high energy cosmic rays and therefore, neutrinos and gamma-rays may be produced in hadronic 
interactions with the surrounding medium. The gamma-ray light curves of blazars measured by 
the LAT instrument on-board the Fermi satellite reveal important time variability information. A 
strong correlation between the gamma-ray and the neutrino fluxes is expected in this scenario. 

An unbinned method based on the minimization of a likelihood ratio was applied to a subsample 
data collected in 2008 (61 days live time). By looking for neutrinos detected in the high state periods 
of the AGN light curve, the sensitivity to these sources has been improved by about a factor 2 with 
respect to a standard time-integrated point source search. First results on the search for ten bright 
and variable Fermi sources are presented. 



I. INTRODUCTION 

Neutrinos are unique messengers to study the high- 
energy universe as they are neutral and stable, in- 
teract weakly and therefore travel directly from their 
point of creation to the Earth without absorption. 
Neutrinos could play an important role in under- 
standing the mechanisms of cosmic ray acceleration 
and their detection from a cosmic source would be 
a direct evidence of the presence of hadronic accel- 
eration. The production of high-energy neutrinos 
has been proposed for several kinds of astrophysical 
sources, such as active galactic nuclei (AGN), gamma- 
ray bursters (GRB), supernova remnants and micro- 
quasars, in which the acceleration of hadrons may oc- 
cur (see Ref. [l| for a review). 

Radio-loud active galactic nuclei with their jets 
pointing almost directly towards the observer, the 
so-called blazars, are particularly attractive potential 
neutrino point sources, since they are among the most 
likely sources of the observed ultra high energy cosmic 
rays and therefore, neutrinos and gamma-rays may be 
produced in hadronic interactions with the surround- 
ing medium. The gamma-ray light curves of blazars 
measured by the LAT instrument on-board the Fermi 
satellite reveal important time variability information 
on timescale of hours to several weeks, with intensi- 
ties always several times larger than the typical flux of 
the source in its quiescent state @|- A strong correla- 
tion between the gamma-ray and the neutrino fluxes 
is expected in this scenario. 

This paper presents the results of the first time- 
dependent search for cosmic neutrino sources by the 
ANTARES telescope. The data sample used in this 



analysis and the comparison to Monte Carlo simula- 
tions are described in Section 2, together with a dis- 
cussion on the systematic uncertainties. The point 
source search algorithm used in this time-dependent 
analysis is explained in Section 3. The search results 
are presented in Section 4 for ten selected candidate 
sources. 



II. ANTARES 

The ANTARES Collaboration completed the con- 
struction of a neutrino telescope in the Mediterranean 
Sea with the connection of its twelfth detector line 
in May 2008 Q. The telescope is located 40 km off 
the Southern coast of France (42°48'N, 6°10'E) at a 
depth of 2475 m. It comprises a three-dimensional ar- 
ray of photomultipliers housed in glass spheres (opti- 
cal modules [4j|), distributed along twelve slender lines 
anchored at the sea bottom and kept taut by a buoy 
at the top. Each line is composed of 25 storeys of 
triplets of optical modules (OMs), each housing one 
10-inch photomultiplier. The lines are subject to the 
sea currents and can change shape and orientation. A 
positioning system based on hydrophones, compasses 
and tiltmeters is used to monitor the detector geome- 
try with an accuracy of 10 cm. The main goal of the 
experiment is to search for high energy neutrinos with 
energies greater than 100 GeV by detecting muons 
produced by the neutrino charged current interaction 
in the vicinity of the detector. Due to the large back- 
ground from downgoing atmospheric muons, the tele- 
scope is optimized for the detection of upgoing muons 
as only they can originate from neutrinos. 

Muons induce the emission of Cherenkov light in 
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the sea water. The arrival time and intensity of the 
Cherenkov Hght on the OMs are digitized into hits and 
transmitted to shore. Events containing muons are 
selected from the continuous deep sea optical back- 
grounds due to natural radioactivity and biolumines- 
cence. A detailed description of the detector and the 
data acquisition is given in Q.The arrival times of 
the hits are calibrated as described in A LI hit 
is defined either as a high-charge hit, or as hits sepa- 
rated by less than 20 ns on OMs of the same storey. At 
least five LI hits are required throughout the detector 
within a time window of 2.2 /xs, with the relative pho- 
ton arrival times being compatible with the light com- 
ing from a relativistic particle. Independently, events 
which have LI hits on two sets of adjacent or next-to- 
adjacent floors are also selected. 

The data used in this analysis were taken in the pe- 
riod from September 6 to December 31, 2008 (54720 
to 54831 modified Julian days, MJD) with the twelve 
line detector. This period overlaps with the availabil- 
ity of the first data from the LAT instrument onboard 
the Fermi satellite. The corresponding effective live 
time is 60.8 days. From the timing and position infor- 
mation of the hits, muon tracks are reconstructed us- 
ing a multi-stage fitting procedure, based on Ref. 7]. 
The initial fitting stages provide the hit selection and 
starting point for the final fit. The final stage con- 
sists of a maximum likelihood fit of the observed hit 
times and includes the contribution of optical back- 
ground hits. Upgoing tracks are also required to have 
a good reconstruction quality. The latter is quanti- 
fied by a parameter, A which is based on the value of 
the likelihood function obtained for the fitted muon 
(see Ref. for details). The angular uncertainty ob- 
tained from the muon track fit is required to be smaller 
than 1 degree. For this analysis, events are selected 
with A > —5.4. This value results in an optimal com- 
promise between the atmospheric neutrino and muon 
background reduction and the efficiency of the cos- 
mic neutrino signal with an assumed spectrum pro- 
portional to E~'^ , where E^, is the neutrino energy, 
which gives the best 5cr discovery potential. The re- 
sulting sample consists of 628 events obtained in 60.8 
days. The simulations indicate that the selected sam- 
ple contains 60% atmospheric neutrinos; the rest be- 
ing mis-reconstructed atmospheric muons. 

The angular resolution of the reconstructed neu- 
trino direction can not be determined directly from 
the data and has to be estimated from simulation. 
However, comparison of data and Monte Carlo in 
which the time accuracy of the hits was degraded by 
up to 3 ns constrains the uncertainty of the angu- 
lar resolution to about 0.1° Q. Figure [T] shows the 
cumulative distribution of the angular difference be- 
tween the reconstructed muon direction and the neu- 
trino direction for an assumed spectrum proportional 
to E~'^. For the considered period, the median reso- 
lution is estimated to be 0.5 ± 0.1 degrees. 




10-^ 10"' 1 10 10^ 



FIG. 1: Cumulative plot of the distribution of the angle 
between the true Monte Carlo neutrino direction and the 
reconstructed muon direction for upgoing neutrino 

events selected for this analysis. 



III. TIME-DEPENDENT SEARCH 
ALGORITHM 

This time dependant point source analysis is per- 
formed using an unbinned method based on a likeli- 
hood ratio maximization. The data is parameterized 
as a mixture of signal and background. The goal is to 
determine, at a given point in the sky and at a given 
time, the relative contribution of each component and 
to calculate the probability to have a signal above a 
given background model. The likelihood ratio A is the 
ratio of the probability density for the hypothesis of 
background and signal (Hsig+bkg) over the probability 
density of only background (H^kg): 

^^^l^ ^Ps^gia^^t^) + (1 - ^)PbkgiS^,t^) 

Pbkg{a^,U) 

(1) 

where n^ig is the unknown number of signal events 
determined by the fit and N is the total number of 
events in the considered data sample. Psigicti,ti) 
and Pbkg{Si,ti) are the probability density functions 
(PDF) for signal and background respectively. For a 
given event i, ti, Si and ai represent the time of the 
event, its declination and the angular separation from 
the source under consideration. 

The probability densities Psig and Pbkg are factor- 
ized into a purely directional and a purely time-related 
component. The shape of the time PDF for the signal 
event is extracted directly from the gamma-ray light 
curve assuming proportionality between the gamma- 
ray and the neutrino fluxes. For signal events, the 
directional PDF is described by the one dimensional 
point spread function (PSF), which is the probabil- 
ity density of reconstructing an event at an angular 
distance a from the true source position. The direc- 
tional and time PDF for the background are derived 
from the data using the observed declination distri- 
bution of the selected events and the observed one- 
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day binned time distribution of all the reconstructed 
muons respectively. Figure [2] shows the time distribu- 
tion of all the reconstructed events and the selected 
upgoing events for this analysis. Once normalized to 
an integral equal to 1, the distribution for all recon- 
structed events is used directly as the time PDF for 
the background. Empty bins in the histograms cor- 
respond to periods with no data taking (i.e. detector 
in maintenance) or with very poor quality data (high 
bioluminescence or bad calibration). 
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FIG. 2: Time distribution of the reconstructed events. 
Upper histogram (black line): distribution of all recon- 
structed events. Bottom filled histogram (red): distribu- 
tion of selected upgoing events. 

The null hypothesis is given with Usig — 0. The ob- 
tained value of Xdata on the data is then compared to 
the distribution of A given the null hypothesis. Large 
values of Xdata compared to the distribution of A for 
the background only reject the null hypothesis with 
a confident level equal to the fraction of the scram- 
bled trials above Xdata- This fraction of trials above 
Xdata is referred as the p-value. The discovery poten- 
tial is then defined as the average number of signal 
events required to achieve a p-value lower than Str in 
50 % of trials. Figure [3] shows the average number of 
events required for a 5cr discovery (50 % C.L.) pro- 
duced in one source located at a declination of -40° 
as a function of the total width of the flare periods. 
These numbers are compared to that obtained with- 
out using the timing information. Using the timing 
information yields to an improvement of the discov- 
ery potential by about a factor 2-3 with respect to a 
standard time- integrated point source search Q. 



IV. SEARCH FOR NEUTRINO EMISSION 
FROM GAMMA-RAY FLARE 

This time-dependent analysis has been applied to 
bright and variable Fermi blazar sources reported 
in the first year Fermi LAT catalogue [l^ and in 
the LB AS catalogue (LAT Bright AGN sample [ll|). 
The sources located in the visible part of the sky 
by Antares from which the averaged 1 day-binned 



FIG. 3: Average number of events (solid line) required 
for a 5(7 discovery (50 % probability) from a single source 
located at a declination of -40° as a function of the width 
of the flare period [at) for the 60.8 day analysis. These 
numbers are compared to that obtained without using the 
timing information (dashed line). 



flux in the high state is greater than 80 10^^ 
photons. cm~^.s~^ above 100 MeV in the studied time 
period and with a significant time variability are se- 
lected. This list includes six flat spectrum radio 
quasars (FSRQ) and four BLlacs. Table U lists the 
characteristics of the ten selected sources. 



Name 


OFGL name 


Class 


RA [°] 


Dec [°] 


Redshift 


PKS0208-512 


J0210.8-5100 


FSRQ 


32.70 


-51.2 


1.003 


AO0235+164 


J0238.6+1636 


BLLac 


39.65 


16.61 


0.940 


PKS0454-234 


J0457. 1-2325 


FSRQ 


74.28 


-23.43 


1.003 


OJ287 


J0855.4-f2009 


BLLac 


133.85 


20.09 


0.306 


WComae 


J1221. 74-28.14 


BLLAc 


185.43 


28.14 


0.102 


3C273 


J1229. 1+0202 


FSRQ 


187.28 


2.05 


0.158 


3C279 


J1256.1-0548 


FSRQ 


194.03 


-5.8 


0.536 


PKS1510-089 


J1512.7-0905 


FSRQ 


228.18 


-9.09 


0.36 


3C454.3 


J2254.0+1609 


FSRQ 


343.50 


16.15 


0.859 


PKS2155-304 


J2158.8-3014 


BLLac 


329.70 


-30.24 


0.116 



TABLE I: List of bright variable Fermi blazars selected for 
this analysis [TTI ]. 



The light curves published on the Fermi web page 
for the monitored sources [l^l are used for this analy- 
sis. They correspond to the one-day binned time evo- 
lution of the average gamma-ray flux above a thresh- 
old of 100 MeV since August 2008. The high state 
periods are defined using a simple and robust method 
based on three main steps. Firstly, the baseline is 
determined with an iterative linear fit. After each 
fit, bins more than two sigma {(Jbl) above the base- 
line (BL) are removed. Secondly, seeds for the high 
state periods are identified by searching for bins sig- 
nificantly above the baseline according to the criteria: 
[F-ap) > {BL+2*aBL)+F > {BL+3*aBL) where F 
and (Tp represent the flux and the uncertainty on this 
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flux for each bin, respectively. For each seed, the adja- 
cent bins for which the emission is compatible with the 
flare are added if they satisfy: {F — ap) > {BL+cjbl)- 
Finally, an additional delay of 0.5 days is added be- 
fore and after the flare in order to take into account 
that the precise time of the flare is not known (1- 
day binned light curve). With this definition, a flare 
has a width of at least two days. With the hypothe- 
sis that the neutrino emission follows the gamma-ray 
emission, the signal time PDF is simply the normal- 
ized light curve of only the high state periods. 

For nine sources, no coincidences are found. For 
3C279, a single high-energy neutrino event is found in 
coincidence. This event is located at 0.56° from the 
source location during a large flare in November 2008. 
The pre-trial p- value is 1.0 %. FigureH]shows the time 
distribution of the Fermi gamma-ray light curve of 
3C279 and the time of the coincident neutrino event. 
This event was reconstructed with 89 hits distributed 
on ten lines with a track fit quality A = —4.4. The 
post-trial probability computed taking into account 
the ten searches is 10 % and is thus compatible with 
background fluctuations. 
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FIG. 4: Gamma-ray light curve (dots) of the blazar 3C279 
measured by the LAT instrument onboard the Fermi satel- 
lite above 100 MeV. The light shaded histogram (blue) in- 
dicated the high state periods. The dashed line (green) 
corresponds to the fitted baseline. The red histogram dis- 
plays the time of the associated ANTARES neutrino event. 



V. SUMMARY 

This paper presents the first time-dependent search 
for cosmic neutrinos using the data taken with the full 
twelve line ANTARES detector during the last four 
months of 2008. For variable sources, time-dependent 
point searches are much more sensitive than time- 
integrated searches for variable sources due to the 
large reduction of the background. This search was 
applied to ten very bright and variable Fermi LAT 
blazars. The most significant correlation was found 
with a flare of 3C279 for which one neutrino event 
was detected in time/direction coincidence with the 
gamma-ray emission. The post-trial probability is 
about 10 %. Upper limits were obtained on the neu- 
trino fiuence for the ten selected sources. 
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